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Abstract. Strongly coupled plasmas in which the interaction energy exceeds the 
kinetic energy play an important role in many astrophysical and laboratory systems 
including compact stars, laser plasmas and dusty plasmas. They exhibit many unusual 
collective properties, such as liquid or crystalline behaviour, peculiar oscillation spectra 
and transport properties. Recently, strongly coupled plasmas were studied in the 
presence of a strong magnetic field by computer simulations, and strong modifications 
of their transport properties and oscillation spectra were observed. While strong 
magnetization is common in stellar systems it is practically impossible to achieve in 
complex plasmas due to the large mass of the dust particles. Here we discuss a recently 
demonstrated approach to achieve very strong "magnetization" by a rotation of the 
neutral gas, and we present new results for macroscopic two-dimensional systems. 
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1. Introduction 

Strongly coupled plasmas in which the Coulomb interaction energy per particle 
exceeds the thermal energy are of high interest in dusty plasmas, trapped ions, 
colloidal systems, quantum confined semiconductor structures as well as in diverse 
astrophysical environments. In recent years much attention has been devoted to study 
the structural properties, oscillation spectra and transport properties of these systems 
both theoretically and experimentally. A particularly detailed analysis was performed 
for colloidal and dusty plasmas, e.g. [1112], where one takes advantage of the large particle 
size (micrometers) that allows for single-particle resolution and optical diagnostics. 
Therefore, the latter plasmas have been successfully used as a unique "laboratory for 
strong correlations" in Coulomb systems in general [2] allowing, for example, to make 
predictions of important structural and dynamical properties of such exotic objects as 
white dwarf stars or neutron stars, which are not directly accessible to experimental 
studies. 

Many strongly correlated plasmas are, in addition, subject to a strong magnetic 
field which may significantly alter the plasma behaviour. Examples are the electron- 
hole plasma in semiconductor quantum wells (quantum Hall systems), magnetic target 
fusion scenarios or the plasma in the outer layers of neutron stars |3]. Recent theoretical 
studies showed that the magnetic field drastically alters the particle and energy transport 
in these systems, in particular the diffusion coefficient see Sec. [2] and the collective 
oscillation spectrum jH E]- It is tempting to directly verify these predictions in the 
laboratory with dusty plasmas, making use of the above mentioned unique spatial and 
temporal resolution. Unfortunately, due to the large particle size, it is practically 
impossible to magnetize a dusty plasma, see Sec. [2j We have recently presented a 
solution to this dilemma: the key idea is to put the dust particles into rotation via a 
rotation of the neutral gas. Then, in the rotating frame, the dust particles are subject 
to the Coriolis force that is equivalent to the Lorentz force of a magnetic field. A brief 
derivation and proof of principle experiment was reported in for a system of 4 dust 
particles in a plane. In this paper we present a more detailed derivation of the equations 
of motion of a rotating dusty plasma, extend the analysis to macroscopic 2D systems and 
analyze the achievable "magnetization" , see Sec. [3] Before that we introduce, in Sec. [2| 
the relevant parameters of a strongly correlated magnetized plasma, discuss important 
special cases and the main effects of the magnetic field. 

2. Magnetization effects in strongly coupled plasmas 

A classical strongly coupled Coulomb system in thermodynamic equilibrium is described 
by a single dimensionless parameter, the coupling parameter 

^ ^ QV47reo 
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Figure 1. Trajectories of a moderately (F = 2, top, five plasma periods) and strongly 
coupled (r — 100, bottom, fifty plasma periods) one-component 2D plasma for different 
magnetic field strengths. Left column: /3 = 0, middle column: (3 = 1, right column: 
/? = 4. Results of MD simulations. 



This assumes that the dynamics can be described by that of an effective one-component 
system (OCP) where the second component assures neutrahzation (and possibly 
screening of the interaction) but does not dynamically participate. This is well fulfilled 
in the case of large mass asymmetry of the components, such as in dusty plasmas 
(if ion streaming is not important) or in dwarf stars and the envelope of neutron 
stars and will be at the basis of the present analysis. Thus, eq. ([T| contains only 
properties of the heavy component: the particle charge Q, the mean inter-particle 
distance (Wigner-Seitz radius) Ows and the temperature T that is related to the thermal 
velocity vt = {kBT/mY^'^, where m is the mass. A second quantity that characterizes 
the strength of the Coulomb interaction is the plasma frequency (dim= 2, 3 is the 
dimensionality) , 

a.p = dimV^./^!/^. (2) 

V "^«ws 

Typical correlation effects are depicted in figure [T| where we show the trajectories of 
several particles computed from a molecular dynamics (MD) simulation. For moderate 
coupling (r = 2, top left) free flight periods are occasionally interrupted by collisions 
that bend the trajectories. In contrast, for strong coupling (F = 100, bottom left), the 
interaction is so strong that particles are trapped in local minima of the total potential 
energy ( "caging" ) exhibiting localized oscillations. 

If a magnetic field is turned on, the particles are, in addition to the Coulomb force, 
subject the Lorentz force giving rise to gyration around the field lines with the cyclotron 
frequency and the Larmor radius 
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Figure 2. Approximate regimes of strongly correlated magnetized plasmas in the 
magnetization-coupling parameter plane. Lines of constant parameter 6 for 
dim = 3] are shown together with important examples of such systems (MTF = 
magnetized target fusion [Slinj). 



The effect on the particle trajectories is seen in figure [T] In the middle column B- 
field and correlations have a comparable effect on the trajectories whereas, in the 
right column, the B-field clearly dominates, giving rise to the familiar helical orbits. 
At strong coupling, evidently, the orbits are much smaller, nevertheless the particles 
remain localized around their mean positions. We may quantify the relative strength 
of correlations and magnetic field by introducing dimensionless energy and length 
parameters. 



/3 = ^ = 5J^^, (4) 
Up V dim m 

flws P V dim ■ 1 

Typical parameter regimes for various strongly coupled plasmas are shown in figure |2] 
From the trajectories it should be expected that macroscopic transport properties 
of the plasma will be strongly influenced both by Coulomb correlations and the magnetic 
field. An increase of F leads to localization and thus tends to reduce the particle mobility. 
Similarly, an increase of the field strength will also tend to confine the particles in the 
direction perpendicular to B. The combined effect of Coulomb correlations and magnetic 
field on the particle transport in a 3D OCP has recently been studied by first-principle 
MD simulations |1] and confirmed these expectations as demonstrated in figure |3} There 
is a systematic reduction of the diffusion coefficient with F, left figure. With increased 
magnetic field the diffusion perpendicular to B is further inhibited (right figure). An 
interesting effect is seen in the longitudinal diffusion coefficient. Although there is no 
component of the Lorentz force parallel to B and particle transport in this direction 
should remain unaffected, this picture breaks down in a strongly correlated plasma. 
The reason is that, upon increase of F, large-angle scattering occurs more frequently 
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which couples the motion perpendicular and parallel to the field. As a result the diffusion 
parallel to B becomes also influenced (and suppressed) by the fleld, and eventually (see 
curves for F = 100) the two diffusion coefficients approach each other. 

We also point out that the decay of both diffusion coefficients approaches, for large 
B-fields, a.l/B power law which is familiar from ideal high-temperature plasmas (Bohm 
diffusion). This indicates that here collective mechanisms such as transport mediated 
by magneto-plasmon excitations become dominant. In fact, the oscillation spectrum 
of magnetized strongly correlated plasmas has recently been computed numerically 
for 2D [5l [To] and 3D systems [6]. The magnetic field was found to drastically alter 
the oscillation spectrum giving rise to nonlinear effects such as Bernstein modes and 
higher harmonics. Again, these effects become relevant for sufficiently strong fields 
characterized by (3 > 0.5. 

Summarizing, we note that correlations and magnetic field (for /3 < 5) suppress 
the particle mobility by about four orders of magnitude compared to the case of an 
unmagnetized weakly coupled plasma which is crucial for a broad range of laboratory 
plasma applications and for many astrophysical problems, such as the evolution of 
stellar objects. It would, therefore, be highly desirable to probe this regime of strongly 
correlated and strongly magnetized plasmas with F ^ 1 and (3 > 0.5 in a laboratory 
experiment with dusty plasmas making use of the unprecedented resolution capabilities, 
e.g. [H E]- After extreme values of F are now routinely achieved, let us now estimate 
the possible values of /3. Assuming a mass of m ^ 10~^^kg, one obtains for a typical 
experiment (ows ~ 2.5 ■ 10""^ m [11]) and a field of B = 4T, common for superconducting 
magnets, a magnetization P ^ 4 ■ 10"'^. This value may be increased if one uses 
smaller (i.e., lighter) particles. However, even then, this is still below the range of 
interest. Furthermore, one would lose the key advantages of dusty plasmas: first, 
smaller particles collect a lower charge which decreases F. Second, the diagnostics 
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of very small particles becomes increasingly difficult, and more sophisticated methods 
other than video microscopy may become necessary to reliably detect the dust grains. 
Finally, while the dust is practically un-magnetized, electrons and (partially) ions are 
magnetized resulting in very different plasma properties. 



3. Rotating dusty plasmas 

A way out of this seemingly hopeless situation has been recently proposed and 
demonstrated in [7]. The key is the fourth component of a complex plasma — the neutral 
gas — which has been completely disregarded so far in our discussion. Even though the 
typical dust-neutral damping rates are very low, in particular compared to colloidal 
suspensions [T], the small yet finite friction force is sufficient to control the global 
motion of the dust component. This was demonstrated in experiments by Carstensen et 
al. |T2], who used a spinning electrode to create a global rotation of the neutral gas in 
the discharge chamber. One application of their method involved the centrifugal force, 
which was used to determine the particle charge [I3]. It was recently demonstrated [7] 
that the second apparent force in a rotating reference frame — the Coriolis force — can be 
used to "magnetize" the dust grains. We now present details of this approach. 



3.1. Effective magnetization 



In the laboratory system, the Langevin equation for the dust grains — being at positions 
Vi = {xi,yi,Zi), interacting via the force Fintifij) = —ViCpi^Pij, Zij) and confined in a 
potential V{p,z) — reads [7j 

N 

mri{t) = -ViV{pi, Zi) + ^ i^mt(n,i) - [ri{t) - u{ri)\ + fi{t). (6) 



Here, z) is the interaction potential, p = a/x^ + y^, Zij = Zi — zj, pij = . . , 
u{r) the flow field of the neutral gas, and u the damping rate. The random force 
fi{t) is characterized by {fi{t)) = and {fi{t)f^{t')) = 2mukBTJij6'''^S{t - 1'), where 
a, /3 E {x,y,z} [neutral gas temperature Tj. For the special case of a uniform angular 
velocity Q, u{r) = (Qcz) x r, one can adopt the point of view of an observer at the 
origin of the coordinate system who co-rotates with the gas column. The coordinates 
in the laboratory {r{t)} and the rotating frame {r{t)} are connected via (e.g., |14j ) 

r(t)=i?(t)f(t), 

r(^) = i^(^)[(^]e,)xf(^) + f(^)], (7) 

f(t) = R{t) [{n e,) X {{n e,) X f{t)} +(2 ^] e,) x f{t) + f{t)\ , 
where the rotation matrix reads 

/ cos(^]^) - sm{Vtt) o\ 
R{t) = sin(fit) cos(fit) . (8) 
\ ij 
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Using ([7]) in ([6]) and multiplying by R^{t) yields [7] 

N 

mf,{t) = -ViV{p„ Zi) + -^mt(nj) + Fcor{f^^) - ymf,{t) + (9) 

where -Fcor(^) = 2m r x (fie^) denotes the Coriolis force. 

Since the angular rotation frequency of the gas is constant, eq. ^ has the same 
form as the Langevin equation for a stationary gas. The confinement potential in the 
rotating frame must be complemented by the centrifugal potential, i.e., 

— Tfl 

v{-prz) = v{-prz)-^nV (10) 

whereas the interaction forces are equivalent to those in the laboratory frame. The 
random force is given by fi{t) = R^(t)fi(t). While = the x and y 

components require more attention. As an example, we consider ff{t) = cos{Qt)f^{t) + 
sm{flt)ff(t) and calculate its properties. Since (fiit)) = {fiif)) = and both /f (t) 
and /f (t) are not correlated with we have {ff{t)) = and {ff{t)fi{t')) = 0, 

respectively. Using {n{t)f^{t')) = and = {f!it)f!it')) = 2mvk^TJ{t - 

t'), one can further show that {ff{t)ff{t')) = and {ff{t)ff{t')) = 2mukBTJ{t - t'). 
The calculations for the y component are analogous, and the results can be summarized 
as {fi{t)) = and {f^{t)f^{t')) = 2mukBTJij6''^6{t - t'), i.e., the mean and the 
correlation function of the random force in the laboratory and the rotating frame are 
equivalent. 

The mathematical equivalence between the Lorentz force -Fl(^) = Qr x [B Bz) and 
the Coriolis force is the reason why the approach can be used to study "magnetized" 
plasmas. When observed in the rotating frame and taking the centrifugal potential into 
account, the dust particles follow the same physical laws as in a constant magnetic field 
Sgff = {2mVL/Q)ez. The important quantity for the magnetization parameter (|4]) is the 
cyclotron frequency Wc, which can now be replaced by 2VL. Since rotation frequencies 
of a few Hz are easily reached, strong magnetization of the dust component becomes 
possible, as was verified in experiments with a small cluster [7]. 

Thus rather than trying to overcome the problems encountered when using 
expensive superconducting magnets, we exploit the advantages a dusty plasma offers. 
The coordinate transformation from the laboratory to the rotating frame is easily 
performed, and the same data analysis techniques as for a stationary plasma can be 
applied. The variation of the plasma parameters (coupling, screening) induced by 
the centrifugal force are negligible for small rotation frequencies and can be directly 
accounted for at higher rotation speeds. While the scaling relations for small systems 
have already been discussed [7j, we now focus on larger ensembles and the transition to 
bulk behaviour. 

3.2. Parameter variation in hulk systems 

Particle confinement in dusty plasma experiments is provided by various forces including 
gravity, electric fields and thermophoresis. In the following, we consider a two- 
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Figure 4. Radial density profile oi N = 1500 particles in a rotating gas flow. Left: 
Radial density for n/LUj_ = 0, 0.7, 0.9 with k(0) = 1 and r(0) = 10. The dashed lines 
show a theoretical result for F — > cx) derived by Totsuji et al. |16j [equation (5.17) in 



their paper], where the screening parameter been adjusted according to (11 1. The 
length scale has been modified accordingly. Right: Comparison of the density profile 
at finite rotation frequencies (lines) with results at ft — (symbols) and adjusted 
screening and coupling parameters. The wiggles at large radii are due to the onset of 
shell formation which is not included in the analytical model |16) . 



dimensional system (implying strong vertical confinement) and assume that the particles 
interact via the Yukawa potential 0(r) = (Q^/47reo r) exp(— r/A), where A denotes the 
Debye length. The effective in-plane confinement is given by a harmonic potential 
V{r) = mu'j_r'^/2, and all higher order terms are neglected. Their contribution may 
become important for large cluster sizes or very high rotation frequencies. 

We first define the dimensionless screening parameter k{Q) = a{Q)/X, the 
coupling parameter T{Q) = Q"^ /{AireQa^Q) k^T^), and the friction coefficient 7(^^) = 
iy/u{n), where a{n) = [QV47reomw2(fi)] V3 ^nd w(fi) = ^ujI - il"^. They are useful 
dimensionless parameters for a finite system and scale as 



V ^l) ' r(o) V ^IJ ' 7(0) V 

The scaling behaviour is illustrated in figure |4| where we show the radial density profile in 
a stationary gas and at finite rotation frequencies obtained from the numerical solution 
of the Langevin equations ^ with = 1500 dust particles [|} As shown in the left 
panel, the density drops as we increase the rotation frequency. If the screening and 



coupling parameters are adjusted according to (11), we find excellent agreement with 
simulations at 1] = (right panel), as expected. The Coriolis force does not affect the 
static properties of the system. 

Experiments typically use large yet finite dust monolayers containing a few 
thousand particles to investigate bulk properties of dusty plasmas [T71 [HI HE]. The 
applicability of this approach was investigated theoretically by Sheridan [19]. To 
give estimates for the scaling of the plasma parameters in bulk systems, we refer to 
the theoretical results obtained by Totsuji et al. [16], who developed a theory for 

X We use the "SLO" integrator derived by Mannella in [TS] 
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Figure 5. Scaling of the coupling parameter, screening parameter, and damping 

The magnetization /3b is shown for (i) 

1500 and (ii) Nk{0) 



constant in a 2D plasma according to ( 12 ) 
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300, 



the average density profile of large clusters. One of their results, using the local 
density approximation and neglecting correlations, is shown in figure |4] and shows good 
agreement with the simulation data. While the theory was derived for F — )■ oo, it still 
provides a good description of the mean density for intermediate coupling. In agreement 
with the results of [16J for crystallized plasmas, it slightly underestimates the density 
in the inner region of the trap. Larger deviations were observed for ^ 1 or k'^ ^ 1. 
In the latter case, the theory could be improved by taking the correlation energy into 
account [16j. 

Based on their theory, Totsuji et al. [16] further give an estimate for the mean 
inter-particle spacing, a„s = A(8fi:^/iV)^/^, where the tc'^'^ scaling was found to be in 
good agreement with simulation results [19]. This enables us to estimate the relevant 
dimensionless parameters for a bulk system ('b') rotating with angular velocity Vl as 
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The exponents are very similar to those in (fTTl), but the variation of the screening and 



coupling parameters with ft is slightly reduced. The trend of eqs. (12) is displayed in 
figure [5] Considering the magnetization Ph{^), the factor N'^^^ indicates that it could 
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be more difficult to achieve strong magnetization in large extended systems than in 
small clusters. For the parameters in figure |4| we obtain rb(0) ~ 37 and ki^{0) ~ 0.27, 
which drop to Fb ~ 21 and increase to Kb ~ 0.48 at Q/u± = 0.95. The magnetization 
at this rotation frequency is /3b ~ 0.45, which is sufficient to observe the magnetic 
field effects discussed in Sec. [2j In particular, this should allow to clearly detect 
the magnetoplasmon and magnetoshear modes [201 |2T] . and could even be sufficient 
to observe higher harmonics of the former |H1 [TU]. The actual scaling behaviour in 
experiments could differ from these predictions due to anharmonic contributions to the 
confinement potential or deviations of the interaction from a Yukawa potential. Further, 
the limitations of the theory must be kept in mind. 

4. Conclusions 

In this paper we have given an overview on the combined effect of strong Coulomb 
correlations and strong magnetization. The main trend is a reduction of particle mobility 
which should play an important role for the transport properties of such plasmas. 
Relevant examples have been presented in figure [2j Dusty plasmas provide a unique 
opportunity to accurately study correlation effects, due the large particle size. At the 
same time, this prevents a magnetization of the dust particles, even in the strongest fields 
superconducting magnets can deliver. Instead, strong "magnetization" of the heavy dust 
particles can be achieved by putting them into rotation by using a rotating neutral gas 
column. The appearing Coriolis force is equivalent to the Lorentz force, and the system 
mimics one in a strong magnetic field that is of the order of 10^ T. This concept was 
recently demonstrated experimentally for a dust cluster of = 4 particles [7j. Here, we 
have performed first-principle Langevin dynamics simulations that confirm the scaling 
of the plasma parameters in a finite system. Based on theoretical predictions for the 
density profile [16] , we have derived the scaling of the relevant plasma parameters with 
the rotation frequency in a bulk system, which should allow for a direct comparison 
with experiments. 
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